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SU(1,1) symmetry of multimode squeezed states 

Z Shaterzadeh-Yazdi, P S Turner and B C Sanders^ 

^Institute for Quantum Information Science, University of Calgar 
2500 University Drive NW, Calgary, Alberta, Canada, T2N IN." 

We show that a class of multimode optical transformations that employ linear optics plus two- 
mode squeezing can be expressed as SU(1,1) operators. These operations are relevant to state- 
of-the-art continuous variable quantum information experiments including quantum state sharing, 
quantum teleportation, and multipartite entangled states. Using this SU(1,1) description of these 
transformations, we obtain a new basis for such transformations that lies in a useful representation 
of this group and lies outside the often-used restriction to Gaussian states. We analyze this basis, 
show its application to a class of transformations, and discuss its extension to more general quantum 
' optical networks. 



PACS numbers: 03.65.Fd, 42.50.Ex, 03.67.-a 



I. INTRODUCTION 



pH , The Heisenberg uncertainty principle ensures that light is noisy at a quantum level. Mathematically, this noise 

■ arises because of non-commutativity of the sine and cosine quadratures of light. In the vacuum state, the noise in 
^ I both quadratures are equal, but the noise can be squeezed - i.e. the noise in one quadrature can be reduced at the 
O^. expense of increasing noise in the complementary quadrature - by propagation through nonlinear optical devices such 

as crystals and gases. In quantum information applications, squeezers, which shift fluctuations from one quadrature to 

■ another [l| , provide an entanglement resource through correlated noise in two modes with these correlations stronger 
J> ' than anything possible in a classical description of light [3, Q . 

, The importance of squeezed light in quantum information is evident in its role in some of the most important 
' experiments in quantum information science: teleportation entanglement swapping tests of local realism and 
' Einstein-Podolsky- Rosen (EPR) paradox , and quantum state sharing 0, Q . These continuous variable quantum 
information experiments involve more than two modes, and entanglement between more than two modes is inherent 
in these experiments. 

Here we analyze multimode squeezed light. Specifically, we are interested in the case of two-mode squeezed light 
(which could be generated by a two-mode squeezer or, alternatively, by two single-mode squeezers in opposite phase 
with their outputs mixed at a beam splitter to yield two-mode squeezed light ^]), with the squeezed field distributed 
over multiple modes via linear optics to share the squeezed light between modes. Mathematically, the transformation 
is described by a matrix, which transforms input state vectors into output state vectors. The choice of basis can 
simplify the mathematical description as well as the calculations themselves. For example, if the inputs are Gaussian 
' states, it is convenient to use a Gaussian basis, and transformations can be fully described in terms of the multivariate 
vector of means and the multivariate covariance matrix ^IQJ . For general input states, the Fock number state basis is 
a useful basis to write general transformations; two alternative representations that transform simply under squeezing 
and linear optics are provided by the Wigner function [ll], and by the position representation 's'l . The problem 
with these approaches for general input states is that the size of the matrix grows linearly in the number of modes. 
On the other hand, the use of group theoretical concepts in quantum optics is becoming more and more prevalent (Tsl . 
fisl . fl6| . Understanding the symmetry of a quantum system and identifying its group properties allow one to use a 
large class of mathematical tools to simplify problems ^17j . Here, we show that an alternative basis can be obtained, 
which relies on determining the symmetries of the network transformation and using the powerful machinery of Lie 
group theory, including the Wigner-Eckart theorem p^ . In our basis, the mathematical description has constant 
size independent of the number of modes. Our approach builds upon the SU(1,1) symmetry inherent in combining 
squeezing with linear optics in two modes. 
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In the absence of squeezing, n-mode linear optical transformations are elements of the compact special unitary 
group SU(n) [191] : with squeezers included in the network, optical network transformations are given by elements 
of the symplectic group Sp(2n,]R) [lol . [20l |. The fact that these transformations are members of the symplectic 
group guarantees the simplicity of describing Gaussian state transformations using just covariance matrices and 
transformations of means |2l| , but this simplicity does not carry over to general non-Gaussian states. 

Here, we show that a broad class of quantum optical networks involving two-mode squeezing and linear optics of 
arbitrarily many modes, can be greatly simplified by exploiting a SU(1,1) =Sp(2,M)c Sp(2n,IR) symmetry in such 
systems. This symmetry allows us to describe the transformation with fixed size independent of the number n of 
modes. All examples mentioned above - quantum teleportation, entanglement swapping, state sharing, and tests of 
local realism - have transformations that are members of this class. We exploit this SU(1,1) symmetry by finding a 
basis that reduces the n-modc Fock states into irreducible representations (irreps) of this group. 

In addition to finding the SU(1,1) symmetry for a broad class of interferometers comprising linear optics and 
two-mode squeezing, and finding a convenient basis in which the transformations have fixed size independent of the 
number n of modes, our work points to new directions in studies of more complex interferometers involving more 
squeezers. The Bloch-Messiah theorem applied to quantum optical networks [13, is a powerful tool to reduce such 
interferometers to all linear optical transformations, followed by single-mode squeezers, followed by all linear optical 
transformations. The approach we establish here may provide an alternative to decomposing such interferometers 
by concatenating interferometers into a larger whole. We discuss this possible future direction for research in the 
Conclusions. 



II. BACKGROUND: TWO-MODE SQUEEZING AND THE SU(1,1) LIE GROUP 



Two-mode squeezed states can be generated either by entangling two independent single-mode squeezed states via 
a 50:50 beamsplitter, or by employing the non-degenerate operation of a non-linear medium in the presence of two 
incoming modes [23 |. The unitary operator describing two-mode squeezing is 



Sabiv) = exp -{{ridb + ifa)h^)/2 



(1) 



with d, h the annihilation operators of the incoming modes and 77 e C the squeezing parameter. This operator gives 
a unitary representation of the SU(1,1) Lie group on the Hilbert space of two modes. As such, it is generated by a 
su(l, 1) Lie algebra given by 



k+ = a^b\ 
which satisfy the commutation relations 

The SU(1,1) Casimir invariant is 

1 



= ab. 



= 2K, 
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(2) 



(3) 



(4) 



Eigenvalues of k^ are used to label the irreps of SU(1,1), and eigenvalues of ATq provide an index for a basis of each 
irrep. Denoting such an orthonormal basis by {\k, /i)}, we have the following SU(1,1) action 



k^\k,y) = k{k-l)\k,y). 



(5) 
(6) 
(7) 



Due to the non-compactness of this group, all unitary irreps are infinite dimensional. There are in fact several different 
series of irreducible representations of SU(1,1) distinguished by the domains of these eigenvalues [l3|. For now, we 
are only interested in the usual positive discrete series of irreps where fc is a non-negative half integer and /i takes 
values fc-|-mforTO = 0,l,2,-- - carried by a Hilbert space denoted by Sl'^ . 

If we label the number of excitations for modes a and b by Ua and n;,, respectively, then k satisfies 



k 



nb\ 



1 



(8) 
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As the relabelling of modes a ^ 6 is physically equivalent to the original labelling, we consider the cases ±(na — Ub) 
to be equivalent irreps (hence the absolute value above). Thus we say that each irrep k > occurs twice, and the 
Hilbert space decomposes as 

Jifa (Si -yfb ^ ® 2^+ © 2^t ® ■ ■ • (9) 

2 2 

Assuming that mode a has more photons than mode b, an arbitrary two-mode state jria, is equivalent to SU(1,1) 
weight states with 

k = ^ , M= ^ • (10) 

As the irrep label k is proportional to the photon number difference, and irrep spaces are invariant under SU(1,1), 
the photon number difference is conserved for two-mode squeezers. In other words, the photon number difference in 
the two modes entering these optical elements equals the photon number difference leaving them. 

Discovering an appropriate realization of a symmetry group enables a clear understanding of a system via the 
mathematical properties that are already known for these groups. We can exploit representation-theoretic machinery, 
such as selection rules and branching rules to facilitate calculations [13, [25l. [26.1 . Output states are then characterized 
by generalized coherent states that are known for groups such as SU(1,1) [27| . 

Often, understanding the symmetry of a system also brings with it much needed physical insight as systems become 
more complicated. However, we will limit the scope of this paper to identifying that symmetry and constructing 
representations for a class of linear optical networks plus squeezers that occur in several important optical quantum 
information protocols as described in the next section. 



III. MOTIVATION: THREE-MODE SQUEEZING AND THE SU(1,1) LIE GROUP 



As mentioned in the Introduction, the motivation for this research stems from recent optical quantum information 
experiments. In Fig. [T] we have included simple schematics of these experiments. As is indicated in the figure, all 
these systems incorporate a two-mode squeezer with one (or both) of the squeezer's output modes mixed with a third 
(and fourth) mode on a 50:50 beam splitter, thereby distributing entanglement across three (or four) modes in the 
network; see Fig. O 

Two-mode squeezing action is given in Eq. ([1]); the beamsplitter action on modes a and b is given by 

I3ab{9, 4>) = exp 9 (a^Se'"^ - aSte"'"^) /2 . (11) 
For the given values of the arguments, the transformation performed by the three-mode component of Fig. [2] is [l^] 

(12) 

for 



^iai'5aibi(2i77) = 5a2ai6i(\/2i?7)^+„^ 



B 



± _ 



and 



2ai&i iy^i-q^ 



exp 



x/2 



ai6i - 0261 



177 
71 



(a\b\ - a\b\ 



(13) 



(14) 



It is not difficult to check that the generators of this transformation satisfy a su(l, 1) algebra with ladder operator 



K+ = ^[a\b\-a\b\ 



(15) 



The specifics of this realization have discussed in our preliminary investigation [28l | . In the next section we generalize 
this result to arbitrarily many modes. 
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FIG. 1: (Colour online) The dashed rectangles show that how three-mode or four-mode squeezed states, which are a key source 
of entanglement, can be generated by distributing a two-mode squeezed state over other modes. Such states are a required 
ingredient for various quantum information tasks such as: (a) quantum state sharing [2^ : (b) quantum teleportation [iTI ]: 
(c) entanglement swapping [||; and (d) tests of local realism and Einstein- Podolosky- Rosen paradox [^. 

IV. RESULTS: MULTIMODE SQUEEZING AND THE SU(1,1) LIE GROUP 



In this section, we show that the three- mode entanglement distributing component given in Fig. [2] can be extended 
to arbitrarily many modes while still being generated by a su(l, 1) algebra. We then analyze the representations of 
this algebra on the multimode Fock space, giving the new basis of SU(1,1) weight states. 

We are considering an optical network (see Fig. that comprises one two-mode squeezer in which one output 
state is mixed via beam splitters between r modes, created by {a|}J^]^, and the other mixed with s modes, created 
t'y ^i\i=\- For simplicity we consider only 50:50 beam splitters with beam splitters for the upper r modes having a 
phase (f> = TT and those for the lower s modes having (j) = — tt. Based on the Baker-Campbell-Hausdorff formula [34|. 
the resulting transformation performed by the multimode network is 



6261 ' 



(16) 



with Sa^Bs being a multimode squeezing operator (see Fig. [3]). This operator can be viewed as a Bogoliubov- 
transformed two-mode squeezer and is given by 



Sa^b^ - cxp [-i(r;i,4 + r?*iJ;st)/2 
where Aj. and are the generahzed boson [pseudo-boson) operators 



E 

1=1 



(-1) 



i-i 



(-1) 



s-l 
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(-1) 



s^ll 



(17) 
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FIG. 2: (Colour online) The 'primitive' three-mode optical network in which tripartite entangled states are produced. Such 
states are the first step towards multipartite entangled states and are applied in various optical quantum schemes some of which 
are shown in Fig. [l] 

given in terms of the original optical modes ai and hi for r, s > 2. Of course, these pseudo-operators satisfy the 
canonical bosonic commutation relations 

= [B„Bl] = i. (19) 

The Hamiltonian generating Sa^b^ is a linear combination of the two operators ArBs and These together 

with their commutators are closed under the commutation relations of a su(l, 1) algebra and thus provide a pseudo- 
two-boson realization of su(l, 1) 

K-^ArB,, k+=AlBl Ko^l{AlAr + BlBs + i). (20) 

The Casimir operator for this realization is 

K^ = \ [{AlAr - BlBsf - i] (21) 

with A\.Ar and BIB^ the pseudo-number operators, Na and Nb, respectively. If we define pseudo-number states such 
that 

NaIua} = nA\nA}, NbItlb} = nslriB}, (22) 

then the entire situation is formally equivalent to the two-mode realization given in the Section ini From Eq. we 
can read off the labels for a basis of SU(1,1) weight states \k, ii) in terms of the pseudo-number states \nAnB}' 

k = 1^{Wa - nB\ + l) , /i = fc + m = i(nA + + 1), (23) 

k = - , m - 0,1,2,. •• , (24) 

riA = "ai + "a2 H ^ "a., ub = rib^ + Uf,^ H h nfc^. (25) 

Therefore, the difference between the total photon number for modes ai to and the total photon number for modes 
bi to bs is conserved. 

All that remains in order to specify the new basis is to give the pseudo-number states in terms of the original modes. 
For {tia}, and analogously for jns}, we obtain 

(A'^V^ 

K} = ^^|0} = ^C-|n), (26) 



where n = (ria-^ , rzaa i ' ' ' j n-a^ ) is a partition of ua- 

r 

^nai=nA, 0<na, eZ, (27) 
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FIG. 3: (Colour online) Equivalence between a typical multimode quantum optical network and a pseudo-two-mode squeezer, 
(a) A multimode optical network comprising a two-mode squeezer (Sat) and several 50:50 beam splitters, (b) A pseudo-two- 
mode squeezer {Sab)- 



and |n) = \na-i,na2, • • ■ ,'^■0^) is a multimode number state. From Eqs. (jl8|) and (|26p we have 



1/2 



n 

Ll=l 



(-1) 



(28) 



It is not difficult to verify that these states are orthonormal, {n/i|r7i^} — SuAmA- 

Take for example the experimentally significant three- mode case depicted in Fig. [21 In this case, the pseudo-number 
states, in terms of the original modes ai and 02, are given by 



^ riA / \ 1/2 

1 ^ (^^l)NA-n.J^A\ 



(29) 



which is normalized because 



— y 

2nA /-^ 



ni=0 



= 1. 



(30) 



Assuming that the number of photons in modes oi and 02 exceeds the photon number in mode 61, the /i"^ three- mode 
SU(1,1) weight state in irrep k, in terms of these three photon number states, is 



, k+n-l 

EM) 



"a, =0 



fc + /i - 1 



1/2 



\k + fl- - l,71ai,^ - A:). 



(31) 
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FIG. 4: (Colour online) A multimode network with more than one two-mode squeezer. Such a network can be built by 
concatenating different parellel sections each being of the type of multimode network considered earlier. The known SU(1,1) 
Clebsch-Gordan and Racah coefficients can be used to couple the corresponding representations and to obtain the resultant 
states. 



However, by this assumption we can only get half of the irreps. To get the other half, one needs to exchange fc + /i — 1 
and 12 — k in the above expression. 



V. DISCUSSION 



This SU(1,1) basis enables us to use a variety of mathematical properties that have already been established for this 
group, thereby significantly facilitating calculations. For instance, the SU(1,1) Clebsch-Gordan coefficients 25, 29, 3§| 
are known for coupling different representations and obtaining the resultant states. These coefficients are very handy 
if one wants to concatenate some of these typical multimode networks 'in parallel' (see Fig. 3]). Furthermore, if the 
concatenation takes place between more than two such optical networks, then the Racah coefficients [3l| give the 
transformations between different orderings of the couplings. The Wigner-Eckart theorem [l^ makes the calculation 
of matrix elements of a tensor operator, or any operator generated from the elements of the algebra ( a linear 
combination thereof, such as the squeezing operator) in this SU(1,1) basis much simpler than a direct approach. 
Generalized coherent states [27], orthogonality, and asymptotic behaviour of the matrix elements are some other 
examples of the well known properties for this group. 

The significance of our result is that it allows for arbitrary input states, in contrast to existing methods which 
usually rely upon Gaussian inputs. This technique can therefore be used to characterize a larger class of output 
states. Moreover, because the algebraic structure is independent of the number of modes, complicated multimode 
linear quantum optical networks with two-mode squeezing, that are of this form, can be dealt with much more easily 
than with, for example, Wigner functions, in which one needs to calculate 2n integrals, or covariance matrices for 
which 4ri^ matrix elements must be calculated, where n is the number of modes. 

These advantages come with a cost, namely a multiplicity of SU(1,1) irreps that grows as the number of modes. For 
r + s modes there are always only two pseudo-boson operators, Ar and Bs, leaving r + s — 2 'directions' in Fock space. 
Finding r + s — 2 orthogonal pseudo-boson operators and using them to label multiplicities would be tedious but not 
impossible for large r or s: this is an open question for further research. Also, we have limited our attention to linear 
quantum optical networks with two-mode squeezing with the specific sequence of beamsplitters given by Eq. p6)) — 
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we considered this form since it is a direct extension of the networks in existing experiments. This SU(1,1) reahzation 
will occur for any Bogoliubov transformation that mixes each output port of the two-mode squeezer separately; that 
is, any Bogoliubov transformation that does not mix a modes and b modes (see Fig.|3]). When these modes are mixed 
'in series' by some multiport optical element a far more complicated situation arises, as exemplified in the next section. 

VI. MORE COMPLICATED SCENARIOS 

Consider the network arising from the reconstruction protocol of the quantum secret sharing experiment in Fig. [T] 
Here, there is an extra two-mode squeezer operating on one mode from the upper set and one mode from the lower 
set of the three-mode squeezing scheme given in Fig. [2j Consequently, we obtain the identity 

Sa2bi {v')Sa2aibi (V2u]) = Sa^a.bi {V, v')Sa2bi {v'), (32) 

where iSaaaibi (?y, ?/) is the exponential of a linear combination of twelve terms. However, by fixing two of the param- 
eters, 

7r 

(j) ^ ^ r]' = s' exp(i(/)') — is', 77 = 77* = 2m7r (33) 

for ?n = 0, 1, 2, • • • , we obtain 

Sa2aibi{v,v') = exp 

This operator is generated by 

kx = -i{aibi ~ a\b\), Ky = -i{d2bi ~ alb\), = i{alai - a2d\), (35) 

which satisfy the commutation relations of su(l, 1). A unitarily equivalent form of these operators has been found 
previously by Sebawe Abdalla et al. [32] . One can unitarily transform these operators using 

U = exp 

in order to diagonalizc in the Fock basis: 

K, = Uk,U^ = - a\ai. (37) 
However, this transformation does not simplify the form of the other operators; in particular, the Casimir operator 

^ UK^ui ^ (at^2 - d\dif - 2(4a2 + d\di + i)b\bi - (0^02 + a|ai + 21) + 2i(a^a|Sf - d2dibl). (38) 

Identifying the appropriate representations of this realization of su(l, 1) carried by the Fock space of three modes is 
quite challenging. For one, the lowering operator is a linear combination of both lowering and raising operators. 
This means that the lowest weight fj, = k states are not Fock states but superpositions of Fock states, a situation that 
does not occur in the standard one- and two-mode realizations. More importantly, can have both positive and 
negative eigenvalues in the Fock basis, which, along with the complicated structure of JK^, suggests that we have left 
the realm of discrete SU(1,1) irreps and must consider the less intuitive continuous irreps that do not have extremal 
weights. The problem of which classes of SU(1,1) irreps are supported on multimode Fock spaces is an interesting one 
deserving further study. Understanding these representations would give us a powerful analytical tool since it would 
give a complete symmetry adapted basis for these and ultimately more complicated optical networks, which would 
aid in the investigation of squeezing and entanglement therein. 

VII. CONCLUSIONS 

We have established a novel basis lying in the discrete representation of SU(1,1). This basis is adapted to physical 
problems of multimode squeezing. Such problems occur in various experimental set-ups for optical quantum infor- 
mation schemes involving a two-mode squeezer and several passive optical elements. This technique facilitates the 
calculation of output states without restricting to 'standard' input states like Gaussians. 



cosh 

.V2 



] (ai&i 



albl) + -^sinh ( '- 



(0,201 



d2d\) 



-^(0261 - 4^1) 



(34) 
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We also showed that some interesting mathematical problems arise in more complicated multimode schemes. These 
problems seem inherent to linear quantum optical networks with squeezers and with more than two mixed modes; e.g. 
it also arises in a four- mode realization of squeezing studied by Bartlett et al. [s^j. This interesting complication opens 
questions about whether, by concatenating optical networks, each with just one two-mode squeezer, the description 
remains within a discrete series representation of SU(1,1). Based on the evidence discussed in previous section, it seems 
unlikely. Such problems are important in considering practical quantum information tasks, in studying squeezing as 
a resource and in similar problems; our work reported here is an important step in this direction. 
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